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Recent measurements performed with some direct dark matter detection experiments, e.g. CDMS- 
II and CoGENT (after DAM A/LIBRA), have unveiled a few events compatible with weakly inter- 
acting massive particles. The preferred mass range is around 10 GeV, with a quite large spin- 



independent cross section of 10 



-10" 



In this paper, we recall that a light dark matter 



particle with dominant couplings to quarks should also generate cosmic-ray antiprotons. Taking 
advantage of recent works constraining the Galactic dark matter mass profile on the one hand and 
on cosmic-ray propagation on the other hand, we point out that considering a thermal annihilation 
cross section for such low mass candidates very likely results in an antiproton flux in tension with 
the current data, which should be taken into account in subsequent studies. 



PACS numbers: 95.35.+d, 95.30.Cq, 96.50.S, 12.60.-i 

The DAMA Collaboration has long claimed the detec- 
tion of an annual modulation in their data [l|, l2| , which, 
if interpreted in terms of dark matter interaction with 
the detector, seems to favor weakly interacting massive 
particles (WIMPs) with light masses (see e.a. Q). More 
recently, the CDMS-II and CoGENT [| Collabora- 
tions, have announced excess events in their data. In- 
terpretations in terms of dark matter were performed in 
e.g. [y] and [7J (see also Q for less conventional mod- 
els), indicating that WIMPs with masses around 10 GeV 
could also explain these measurements. In Ref. pj, it was 
notably shown that some of the favored regions may not 
be compatible among each other. 

Interesting constraints on such light WIMPs may ac- 
tually come from colliders [9h11J. From the astrophysi- 
cal point of view, the annihilation of such light WIMPs 
may generate high gamma-ray fluxes which are at the 
edge of exclusion with current measurements [12| . but 
there are still large uncertainties coming from our incom- 
plete knowledge of the detailed dark matter distribution, 
in particular in the centers of galaxies. On the Galac- 
tic scale, cosmic-ray antiprotons also provide interesting 
constraints [1 31 ] since predictions arc less sensitive to the 
choice of the halo profile, but more to the local density 
(see e.g. [14[) — indeed, the relevant annihilation yield is 
averaged over a volume set by the diffusion scale, which 
of the order of a few kpc about the Earth. It was already 
noticed in |15[ that light neutralinos with masses below 
10 GeV might generate antiproton fluxes overshooting 
the data if the dominant annihilation proceeds into bb, 
because the l/m x 2 flux suppression is no longer efficient 
with respect to heavier WIMPs. In their analysis, these 
authors used a smooth cored isothermal halo profile for 
the dark matter distribution in the Galaxy, which was 
not meant to be in clean agreement with the Galactic 
rotation curves. Nevertheless, the largest uncertainties 
did actually come from propagation, notably from the 
relative freedom in setting the vertical extent L of the 
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diffusion zone. Indeed, in their minimal case of L = 1 
kpc, at the pessimistic edge of cosmic-ray nuclei con- 
straints (la . Il7| , the dark matter contribution was shown 
to be dramatically decreased by almost 1 order of mag- 
nitude with respect to the best-fit propagation setup. 

Here we take advantage of the recent works performed 
(i) by Catena and Ullio [l9| (CU10 hereafter) on con- 
straining the Galactic dark matter distribution from 
kinematic data on the one hand, and (ii) by Putzc, 
Derome, and Maurin |20| (PDM10 hereafter) on cosmic- 
ray propagation on the other hand, to improve the an- 
tiproton analysis. 

CU10 notably showed that one could reach interest- 
ing constraints on the local dark matter density by us- 
ing updated tracers of the Galactic dynamics, provided 
some initial assumptions about the dark matter pro- 
file. These assumptions can be made on well-motivated 
theoretical grounds, since the highest-resolution cosmo- 
logical N-body simulations to date focused on Milky- 
Way-like galaxies now seem to converge towards sim- 
ilar predictions, in between an Einasto |2l], [22( and 
an Navarro-Frenk- White [H (NFW) profile (e.g. 0- 
[26j). Using these assumptions, CU10 derived a lo- 
cal dark matter density of p(r Q — 8.25 ± 0.29 kpc) = 
0.386 ± 0.027 GeV/cm 3 in the former case and of p(r = 
8.28 ± 0.29 kpc) = 0.389 ± 0.025 GcV/cm 3 in the lat- 
ter case. We note that these results for the local den- 
sity were confirmed independently by [271 ] . with a com- 
pletely different method — these authors derived p© = 
0.41 ±0.11 GeV/cm 3 . There are still, obviously, large un- 
certainties with respect to the dark matter distribution 
in the inner kpc about the Galactic center, since either 
the baryons and the central black-hole may play impor- 
tant roles there, increasing or decreasing the inner den- 
sity depending on the hypotheses [28|, [29|. Nevertheless, 
we underline that decreasing the density in the inner kpc 
has much less impact on the antiproton flux predictions 
than on the gamma-ray flux predictions because of spa- 
tial diffusion, as it is demonstrated further below. This 
is strengthened by the fact that the (anti)proton propa- 
gation scale increases with energy: antiprotons are more 



UJ10> 



halo model 

CU10a: NFW 
CU10b: Einasto 



I ' r 

n^bb(IOGeV) 




BESS 95-97 
» BESS 98 

* BESS 99 

* BESS 00 
a BESS 02 
a AMS 98 

* PAMELA 10 
— background 



halo model 

CU10a:NFW 
CU105: NFW+core 




kinetic energy E [GeV] 



kinetic energy E [GeV] 



FIG. 1. Left: Predictions of the primary antiproton flux for a 10 GeV WIMP annihilating into bb pairs, for the NFW (CUlOa) 
and the Einasto (CUlOb) profiles derived in [ljj]. The solid black curve is the secondary background predicted using the same 
propagation setup [13]. Right: Impact (i) of imposing a core to the NFW case (CUlOa versus CUlOa, dashed curve), (ii) of 
increasing the solar modulation force field (dotted-dashed curves), and (iii) of decreasing the vertical halo boundary L from 4 
to 3 kpc (dotted curve). 



local if belonging to the energy range of interest here, say 
below 10 GeV, than those of higher energies. Typically, 
the propagation scale A is set by the ratio of diffusion to 
convection and spallation, and is of order of a very few 
kpc for a 1 GeV antiproton, whatever X < L [30] • Other- 
wise, L provides an extra-limit to the propagation scale: 
the probability to escape the diffusion halo strongly in- 
creases for A > L, so that antiprotons can hardly come 
from regions distant by more than a very few times L. 

Concerning uncertainties in the cosmic-ray propaga- 
tion parameters, L, as we stressed above, has the most 
dramatic impact on dark matter signals. Nevertheless, 
the recent MCMC analysis performed by PDM10 showed 
that adopting a very small diffusive halo with L down to 
1 kpc makes it very difficult to accommodate the cur- 
rent cosmic-ray nuclei data. Moreover, large diffusion 
halo models, as large as ~ 10 kpc, are also preferred in 
astrophysical studies of the high-latitude diffu se g amma- 
rav emission measured by Fermi 13 ill (see e.g. [321]). This 
minimal L of 1 kpc proposed in [l4| was actually not re- 
ally motivated by observational facts, but by the sake of 
a very conservative approach. As independent hints in 
this respect, it is indeed commonly observed that the ra- 
dio halos of nearby spiral galaxies have sizes larger than 
1 kpc (e.g. [33|), and a few kpc vertical extent is also fa- 
vored by studies of the Galactic magnetic field (e.g. [34|). 

In the following, we adopt the best-fit model derived 
in [111, which is astonishingly close to the best- fit setup 
derived in PDM10, in which the diffusion halo has a ver- 
tical extent of L = 4 kpc, still much lower than what is 



suggested by the diffuse gamma-ray interpretation. We 
study a 10 GeV dark matter particle candidate, Majo- 
rana fcrmion or scalar, with a thermal annihilation cross 
section (av) =3x 10 _26 cm 3 /s, entirely annihilating into 
bb quark pairs — using other quark flavors would barely 
change the antiproton production; furthermore, taking 
a lower branching ratio into quarks would translate lin- 
early into our analysis results. We have derived the injec- 
tion antiproton spectrum with the public code PYTHIA 
[35| . For the dark matter profile, we consider the NFW 
and Einasto profiles constrained in CU10 and discussed 
above, which we respectively denote CUlOa and CUlOb 
in the following. We also investigate the potential effect 
of imposing a 1 kpc core to the NFW case, which we 
denote CUlOa. 

In Fig. [TJ we show the antiproton flux predictions at 
the Earth obtained with the series of ingredients intro- 
duced above. We also plot the secondary antiproton 
background consistently derived within the same prop- 
agation setup [18[. The data points are taken from [36r - 
|41| . In the left panel, it clearly appears that the differ- 
ence coming from using different halo profiles is negligible 
- at the order of a few percents, hard to see from the 
plot. This comes from the fact that the local normaliza- 
tion of the dark matter density is quite the same, and 
that the global shape does not differ significantly on the 
kpc scale around the Earth. We note that the primary 
contribution originating from dark matter annihilation 
does exceed the secondary background^ below ~ 2 GeV, 
overshooting it by a factor up to almost 5 around 200 



MeV. In the right panel, we illustrate the effect of im- 
posing a core to the NFW case (CUlOa, dashed curve), 
and demonstrate that this has a very poor impact. This 
is due to the limited propagation scale that characterizes 
the transport of low-energy antiprotons. We also study 
the influence of modifying the force field applied for so- 
lar modulation [42[ , increasing it from <fi = 600 MV (solid 
curves) to 900 MV (dotted curves) . Such a change applies 
to both the signal and the secondary background, which 
makes the argument still valid in the strong solar activity 
regime. Finally, to allow a more conservative view, we 
consider a decrease of L down to 3 kpc self-consistently 
with the cosmic-ray nuclei constraints [18|] (dotted curve) 
— such a change has no effect at all on the background 
prediction. We see that even in that case, the predicted 
primary flux exceeds the secondary background by a fac- 
tor of 2, leading again to serious tensions with the data. 
We still emphasize that many independent hints favor a 
large diffusion halo model with L > 3 kpc, as already 
mentioned above. 

Therefore, a light dark matter particle in the ~ 10 GeV 
mass range, annihilating into quark pairs at the thermal 
rate set by the relic density constraints, is expected to 
harden the antiproton spectrum below a few GeV. This 
very likely leads to important tensions with the current 
data. Further accounting for the predicted presence of 
subhalos would make this statement even slightly more 
severe [43| . as well as considering a dark matter disk 
due to subhalo tidal streams trapped into the Galactic 



disk [44H46I ] . Turning the argument around, observing 
a net inflexion in the antiproton spectrum around a few 
GeV could be a hint pointing towards the contribution 
of light WIMPs — it does not seem to be the case in the 
available data. A loophole is of course possible so as to 
decrease the primary signal and escape these constraints: 
e.g. combining cither a lower branching ratio into quarks 
or a lower annihilation cross section with playing with the 
astrophysical parameters, or, for complete safety, taking 
a WIMP mass less than the (anti)proton mass |47| . Con- 
cerning the astrophysical parameters, the ongoing efforts 
to constrain them have started to allow for less freedom 
in the predictions, at least in the domain of local charged 
cosmic rays. 

To conclude, we emphasize that light dark matter can- 
didates considered in the interpretation of direct detec- 
tion signals should be checked against the cosmic-ray an- 
tiproton data, at least whenever their couplings to quarks 
are significant. A more systematic study of this comple- 
mentarity is on-going for different particle physics sce- 
narios beyond the standard model and is about to show 
that some models are already excluded, except in con- 
trived astrophysical situations [48J . 
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